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VARYING  BOUNDARY  CONDITIONS 
WITH  LARGE  DIFFUSIVITY 

by 

Jack  K.  Hale  and  Carlos  Rocha 

ABSTRACT 

For  systems  of  semi linear  parabolic  partial  differential  equations 
on  bounded  domains  with  large  diffusivity  and  homogeneous  boundary 
conditions  close  to  the  Neumann  conditions,  we  associate  a  system  of 
ordinary  differential  equations  (ode's)  from  which  the  dynamics  of  the 
original  system  can  be  inferred.  Small  perturbations  of  the  Neumann 
case  produce  large  perturbations  in  the  ode's  with  corresponding  effects 
on  the  dynamics  of  the  system.  The  same  theory  is  valid  for  functional 
differential  equations.  Applications  are  considered  in  models  for 
control  by  genetic  repression  of  biological  material  in  cells. 


1.  Introduction. 


Consider  the  system  of  parabolic  partial  differential 
equations  (PDE) 

(1.1)  3u/9t  *  Dau  +  f(u);  X  6  fi 

(1.2)  3u/3n  =  E(x)u;  x  e  3fi 

N  n 

where  u  €  R  ,  c  F  ,  n  <  3,  is  a  bounded  open  set  with  3fi 
smooth,  D  =  diag(di ,. . .  .dj^),  E  =  diag(ei ,. . .  .Oj^)  where  each 
d.  >  0  is  constant  and  each  e-:  3n  ->■  F  is  continuous.  Also, 

J  w 

N  N  11 

suppose  that  f:  F  ^  F  is  a  C  ’  -function;  that  is,  is  con¬ 
tinuous  and  has  a  Lipschitz  continuous  first  derivative. 

Our  objective  is  to  study  the  behavior  of  solutions  of  (1.1), 

(1.2)  when  D"\  E  are  small;  that  is,  the  boundary  conditions 
are  close  to  homogeneous  Neumann  conditions  and  the  diffusivity 
is  large. 

2  N 

To  state  the  results,  we  need  some  terminology.  If  X  =  L  (n,F  ), 

D(A)  =  {(;  €  W^*^(n,F^):  3(i>/3n  =  E;  on  3n} 

then  A^)  =  -Da6  is  a  sectorial  operator  and  one  can  define  the 
fractional  powers  A“  of  A,  0  <  a  and  the  space  =  D(A°')  with 
the  graph  norm.  If  3/4  <  a  <  1,  then  X'^  c  W^’^(n,FN)  n  L“(r.,F^) 
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with  continuous  inclusion.  One  can  then  show  that  (1.1),  (1.2) 
defines  a  local  C’  -semigroup  Tp  ^(t)  on  X  for  3/4  <  a  <  1 
(see,  for  example,  Henry  [9,p.75]). 

For  any  set  B  c  X“,  the  u-limit  set  u(B)  of  B  is  defined 
as  a>(B)  *  C£  £(t)B.  A  set  Be  x“  is  said  to  be  in¬ 

variant  if,  for  any  <(i  €  B,  one  can  define  £(t)(t  for  t  €  F 
and  Tp  ^(t)$  €  B  for  t  €  F.  A  set  A  c  x“  is  said  to  be  a 
compact  attractor  for  (1.1),  (1.2)  if  A  is  compact,  invariant  and 
there  is  a  neighborhood  B  of  A  such  that  u.'(B)  c  A. 

Under  the  assumption  that  the  ordinary  differential  equation 

(1.3)  du/dt  =  f(u) 

has  a  compact  attractor  A,  it  was  shown  in  Hale  [5]  that  A  is  a 
compact  attractor  for  (1.1),  (1.2)  if  E*0  provided  that  d>  is 
sufficiently  large,  where  d  »  min(d^,. . .  ,dj^)  and  X  is  the  first 
eigenvalue  of  -l  with  homogeneous  Neumann  conditions.  In  other  words, 
(1.1),  (1.2)  with  E  =  0  behaves  qualitatively  as  the  ODE  (1.3)  if 
dx  is  large. 

We  will  obtain  an  appropriate  generalization  of  this  result  when 
D“\  E  are  small;  that  is,  the  qualitative  properties  of  the  flow  for 
(1.1),  (1.2)  can  be  determined  from  an  ODE  under  certain  hypotheses. 

At  first,  it  is  instructive  to  guess  the  appropriate  ODE.  If 
i7(t)  *  if.  |',u(t,x)dx  and  u(t,x)  is  a  solution  of  (1.1),  (1.2), 
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then,  by  integrating  (1.1)  and  using  the  boundary  conditions  (1.2), 
one  obtains 


du  _ 
3t  ■ 


D|n 


^  E(x)u(t,x)dx  +  |£2l 


-1 


* 

f(u(t,x))dx. 

n 


If  we  assume  that 


E(x)  =  E(x.D). 

(1-4)  D|fiI"^E(x)  -»•  E(x),  E(x)  continuous, 

* 

E(x)dx  =  ? 

and  u(t,x)  -  u(t)  0  as  d  =  min(dp. . .  .d^^)  -*•  »,  then  v  =  lim^^ 
u  satisfies  the  (ODE) 

('■5)  tv  +  f(v)  '*6*' g(v,i;) 

The  "perturbation"  E(x,D)  in  the  boundary  conditions  leads  to  the 
perturbation  in  the  vector  field  f(v)  in  (1.3). 

Under  the  assumption  (1.4),  there  are  N  eigenvalues  x^(D), 
...,X|^(D)  of  Da  with  boundary  conditions  (1.2)  and  corresponding 
unit  eigenfunctions  (D) ,. . .  ,(tj^(D)  such  that 

Ad  “  di ag ( Xi  ( D) , . . . , Xj^( D) )  ■+  ^ 

”  {  =  -| (D)>- •  •  »<;nj(D))  I  =  identity 


as  d  Decompose  as 

X®  *  'J  ® 

a 

U  *  {*  *  ♦jjv:  V  €  R^} 

i/  =  {il,  €  X®:  <i^.<?>*=0  for  all  41  €  u} 

T 

where  <i|/,4»  *  il  (x)4(x)dx. 

h 

If  u(t,x)  is  a  solution  of  (1.1),  (1.2)  and 


(1.6) 

u(t,x)  =  +  w(t,x). 

w(t,*)  € 

a 

then 

(1.7a) 

dv/dt  *  Ap  V  + 

•'n 

f(«P  v+w) 

(1.7b) 

3w/3t  =  Daw  +  f(«p  v+w)  -  4.^1 

$p  f(fp  v+w)  on 

3w/3n  =  Ew  on  dn 


The  main  result  of  the  paper  is  the  following. 

Theorem  1.1. 

Suppose  D,E(*,D)  satisfy  (1.4)  and  the  ODE  (1.5)  has  a  conpact 

N 

attractor  A.  ,  with  V  c  F  being  a  neighborhood  of  A,  , 


such  thet 
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u)(V)  c  A  Then,  for  any  neighborhood  W  of  A  ,  with  C£W  c  V, 

S  »  »  C  »  * 

there  is  a  constant  >  0  and  a  -function  h:  W  {D:  d  >  dp) 

such  that  h(v,D)  -*•0  a^  d  -►  »  and  the  set 

(1.8)  Mp  =  {u  =  4>q  V  +  h(v,D),  V  €  W} 


is  an  exponentially  asymptotically  stable  invariant  manifold  for 
(1.1),  (1.2).  Furthennore,  the  flow  on  this  manifold  is  given  by 
u(t,x)  =  $p(x)v(t)  +  h(v(t),D)(x)  where  v(t)  is  a  solution  of 
the  ODE 


(1.9)  dv/dt  =  Aq  V  + 


f($p  V  +  h(v,D)) 


Apv 


+  G(D,v) 


The  proof  of  this  theorem  is  given  in  Section  2.  It  is 
interesting  to  note  some  implications  of  this  result. 

Since  Ap  ^P  ->  lnj”'^I,h(v,0)  -<•  0  as  d  ->■  relation  (1.6) 
shows  that  the  average  iF  of  u  approaches  jnf'^v  and  v  satisfies 
the  ode 

(1.10)  dv/dt  =  cv  +  ln|  ^fdfil’S) 

A  rescaling  of  v  -*■  ifij^'^v  yields  (1.5). 

The  compact  attractor  A  ^  is  upper  semi  continuous  in  the 
Hausdorff  metric  with  respect  to  -perturbations  of  the  vector 


field  +  f(v)  (see,  for  example.  Hale  [6]).  The  flow  defined 
by  (1.9)  is  equivalent  to  the  flow  defined  by  rescaling  v 
The  new  vector  field  AqV  +  |fi]**G(D, approaches  the  vector 
field  g(v,c)  =  ?v  +  f(v)  in  the  topology  in  W  as  d  -►  «. 

Thus,  (1.9)  has  a  compact  attractor  f  W  and 
cA  f  as  d-^®.  If  the  flow  of  (1.5)  on  A  ^  is  structurally 
stable  in  W,  then  lim^  ^^A^  ^  f  “  f  defined  by 

(1.9)  is  equivalent  to  the  flow  defined  by  (1.5)  in  W.  This  is 
summarized  in 


Corollary  1.2.  For  d  >  dg,  equation  (1.9)  has  a  compact  attractor 
A,,  ^  in  W,  limj  Ar^  _  ^  c  A_  If  the  flow  defined  by  (1.5) 

D,?,t  —  a-*®  u,;,t  ^ 

is  structurally  stable  on  A.  ^  in  W,  then  the  flow  defined  b^ 
(1,9)  is  structurally  stable  on  Ag  ^  f  lH  W  and  is  equivalent 
^  (1 .5)  i£  W. 

The  following  observation  is  also  important.  Although  the 
perturbation  E(-,D)  -*■  0  as  d  ->■  ®,  the  behavior  of  the  flow  de¬ 
fined  by  (1.1)  cannot  be  considered  as  a  small  perturbation  of  the 
flow  defined  by  Neumann  boundary  conditions.  In  fact,  the  matrix 
;  can  be  very  large  and  therefore  the  flows  defined  by  (1.3)  and 
(1.5)  can  be  very  different.  An  illustration  will  be  given  in 
Section  3.  The  flow  defined  by  (1.1)  is  close  to  the  flow  defined 
by  the  equation 
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=  Dau  +  f(u),  X  e  fi 

1^  -  lafil'"'  cu  =  0,  X  e 

for  d  large.  The  similarity  with  Neumann  boundary  conditions  is 
that  the  eigenfunctions  corresponding  to  the  eigenvalues  close  to 
c  are  aporoximately  constant  functions  for  d  large  and  so  the 
average  of  u  satisfies  (1.5). 

Equation  (1.1)  may  also  contain  other  physical  parameters  which 

vary  over  some  set.  In  such  a  case,  equation  (1.5)  would  also  contain 

these  parameters  and  the  flow  may  undergo  various  types  of  bifurcations 

as  the  parameters  vary.  It  would  then  be  of  interest  to  relate  these 

bifurcations  for  the  limit  equation  to  bifurcation  of  (1.9)  for  the 

flow  of  (1.1)  on  the  attractor.  For  example,  suppose  f  =  f(u,r), 

r  e  F,  and  the  limit  equation  (1.5)  with  f  replaced  by  f(u,r)  has 

an  attractor  A  which  is  a  point  v^  for  r  <  0  and  there  is  a 

generic  supercritical  Hopf  bifurcation  at  r  =  0;  that  is,  the  attractor 

A  is  a  disk  for  r  >  0  with  the  boundary  of  A  being  a  stable 
^  ,r  ' 

periodic  orbit  which  attracts  all  points  on  A^  ^  except  the  equilibrium 
point.  Equation  (1.9)  will  also  contain  r, 

(1.9)^  dv/dt  =  ApV  +  G(D,v,r) 

where  ,  G(D,v,r)  ^  f(v,r)  as  d ®  uniformly  for  ir:<rQ 

and  V  e  V  for  some  fixed  >  0  and  neighborhood  V  of  'J;  y«' ,r’ 


This  convergence  will  be  in  the  C  topology  if  f  is  C  *  . 


Furthermore,  the  attractor  „  satisfies  lim.  ^  A 

u,c,r  a-**  D,;,r  ^,r 

1*^1  1  *^0'  remarks  are  immediate  consequences  of  the  proof  of 

Theorem  1.1  below.  Therefore,  if  we  fix  d  i  dQ,  with  dQ  sufficiently 
large,  the  one  parameter  family  of  flows  defined  by  the  vector  fields 
in  (1.9)^,  lr|  1  rp,  have  the  property  that  there  is  an  f(D),  |f(D)l 
<  rQ,  f(D)  0  as  d  -►  »,  such  that  the  attractor  Ap  ^  ^  is  a  point 
for  -rQ  ±r  ^  f(D)  and  is  a  disk  for  f(D)  <  r  £  rp  with  a  generic 
Hopf  bifurcation  occuring  at  f  =  r(D).  The  same  remark  applies  to 
other  elementary  bifurcations  that  occur  in  one  parameter  families. 

It  will  be  apparent  below  that  the  proof  of  Theorem  1.1  is  easily 
adapted  to  the  following  situation.  Let  A  =  A(x)  be  a  positive 
definite,  continuous  nxn  matrix  function  uniformly  for  x  €  n  and 
let  u^  be  the  transpose  of  an  N-dimensional  column  vector.  Consider 
the  equation 


TnT 


(1.10) 


u^  =  D[div  A7u  ]  +  f(x,u)  in 


[Avu^]^n  =  E(u)  in  an 


11  N  N 

where  n  is  the  outward  unit  normal  vector  to  an,  f  €  c  ’  (n  x  p”,  p’) 
and 


E(u) 


e(x,y)u(y)dy 


e(x,y)  *  diag(e.|(x,y),...,enj(x,y)) 


(1.11) 


where  e,  e  are  continuous  on  dn  x  3(2.  Let 


(1.12)  F(u)  =  f  f(x,u)dx. 

■’(2 

Then  Theorem  1.1  Is  valid  with  (1.4)  replaced  by  (1.11),  (1.1), 
(1.2)  replaced  by  (1.10)  and  (1.5)  replaced  by 

(1.13)  u  =  ^u  +  F(u) 

It  will  be  clear  also  form  the  proofs  given  below  that  Theorem  1.1 
remains  valid  even  if  the  equation  (1.1),  (1.2)  is  a  functional  dif¬ 
ferential  equation;  that  is,  f:  C([-r,0],F*'^)  -<•  where  r  >  0  is 
a  given  constant.  Also,  one  can  consider  the  equation  (1.1),  (1.2) 
coupled  with  an  ordinary  differential  equation,  and  obtain  the  version 
of  (1.5)  coupled  with  the  same  ordinary  differential  equation.  This 
is  illustrated  in  Section  4  for  a  problem  in  genetics  of  Busenberg 
and  Mahaffy,  [1].  One  could  also  have  a  functional  differential 
equation  in  (1.10). 


L*j 


2.  Proof  of  Theorem  1.1 


We  will  use  the  notation  Introduced  in  Section  1.  For  a 
fixed  constant  b  >  0  and  a  given  continuous  function  a:  dn  F, 
consider  the  eigenvalue  problem 


(2.1) 


in  n 


3i/3n  =  -a(x)i  on  3n. 


Suppose  a(x)  =  a(x,b)  and 


(2.2) 


b  a(x,b)  Bi(x)  as  b  « 

6l(x)  continuous  B,  = 

'  J  so  ' 


The  smallest  eigenvalue  x.|(b,a)  of  this  problem  is  characterized 
by  the  variational  principle  (see  Courant  and  Hilbert  [2  ,p.  398]) 

(2.3)  -x,(b,a)  =  min{  b!v<j|^dx  +  b  ai^dx:  (^^dx  =  1} 

•’n  ■’3r2  ■'0. 


Also,  if  e.j(b,a)  minimizes  this  functional,  then  c.|(b,a)  is  an 
eigenfunction  corresponding  to  x.|(b,a). 

Lemma  2.1.  If  (2.3)  is  satisfied,  and  )..(b),  k  =  1,2,...  are  the 


eigenvalues  of  (2.1),  then  there  are  constants  >  0,  v  >  0  such  that, 


-n- 

X|^(b)  -yb,  k  ^2,  for  b  >  bg. 

Also,  there  is  a  normalized  eigenfunction  (j)^(b)  such  that 
4ii(b)  -►  l^r**  as_  b  -*•  ». 

Proof:  For  any  e  >  0  there  is  a  bp  >  0  such  that  b  >  bp  implies 

sup  |ba(x)  -  6,(x)l  <  e 

X 

Taking  u(x)  =  it  follows  from  (2.2)  that 

-X,  <  ba  *  Inl'M  (ea-B,)  +  o  <  cj+ek 

hn  hn  ' 

where  K  =  In]"^  |3f2l . 

The  trace  theorem  for  W^*^(n)  implies  there  is  a  constant  Cq 
such  that 

lCQ[f  \V^\^  +  4,^]. 

hQ  ^  Jfi  JQ 

If  we  let  *  sun^|B.|(x)l  and  if.|(b,x)  be  an  eigenfunction  for 
Xi{b),  <t?(b,«)  =  1 ,  then 


Suppose  now  there  is  a  sequence  b.  -►  »  as  j  »  such  that 

J 

wgtbj)  -»■  0  as  j  If  <)>2{bj,*)  is  the  corresponding  eigen¬ 

function, 
and  Inl”^ 


=  1,  then,  as  above,  one  concludes  that  |V(j;2l  0 

1  ^  hi 

<t>2  ^  3S  b  -*■  -<«>.  Thus,  4>2{b,*)  -*•  0  as  b  ®  which 


contradicts  the  fact  that 


=  1.  Thus,  there  is  a  bf.  >  0  and 


a  u  >  0  such  that  ^2^^)  Iv  for  b  ^  bQ.  This  implies  that 
-X2(b)  >  ub  for  b  >  bg.  This  completes  the  proof  of  the  lemma, 


Now  let  us  consider  the  parabolic  PDE 


(2.4) 


3u/8t  =  bAu  in  n 

3u/3n  =  -a(x,b)u  in  sn 


Let  Y  =  L^(n,F),0(A)  =  €  W^’^{n,R):  3(i/3n  «  -a(x,b)c  on  sn} 

and  A(}i  *  -bA(t,  <()  €  D(A).  Then 

(Au)(x)v(x)dx  =  b  [  vu(x)*vv(x)dx  +  b  [  a(x,b)u{x)v(x)dx 

h  'n  hn 

for  u  €  D(A),  V  €  W^’^(fi,F)  and  A  is  formally  self-adjoint. 

1  2 

It  has  a  self-adjoint  extension  to  W  ’  (n,R)  .  Furthermore, 
o(-A)  consists  of  the  eigenvalues  of  (2.1). 

Let  Y®,  0  _<  a  <  1,  denote  the  Banach  spaces  associated  with 
the  fractional  powers  of  A.  Then  (2.4)  generates  an  analytic  semi¬ 
group  e’^^  in  Y“.  For  t  >  0, 
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where  r  is  a  contour  in  the  resolvent  set  p(-A)  with  arg  x 

-*•  ±e  as  |xl  •  for  some  constant  e  e  (-n/Zt-n). 

Let  x^(b)  be  the  first  eigenvalue  of  (2.1)  and  let  (|)i(x) 

»  (^iCb.x)  be  the  corresponding  normalized  eigenfunction-  Let 
'  fl 

<({,,0>  =  (^(x)iii{x)dx  and  let 

Jo 

Y“  *  Y  ^  Y-*-,  Y  *  sp[(t,3. 

a  I 

If  u  e  Y’*,  let  u  =  v^^i  +  w,  V  =  <i;i,u>.  If  u  is  a  solution  of 

(2.4) ,  then 

dv/dt  *  Xi(b)v 
3w/3t  =  bAw  in  n 

and  w  satisfies  the  boundary  conditions 

(2.5)  3w/3t  =  -a(x,b)w  on  3r^ 

Let  us  define  A  =  b"^A|Y^  as  Aw  =  aw  with  the  boundary  conditions 

(2.5)  for  w  e  D(A)  =  D(A)  n  Y'*'.  Then  A  is  sectorial  in  Y'*’,  and, 
by  Lemma  2.1,  c(A)  >  v  Using  [9  ,Theorem  1.5.4]  we  have  for  w  €  Y;*" 
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Ik  t'V"*|w|^2.  t>0 
l*'**w|  „  i  k  e'^^lwl  ,  t  >  0 

yCl  yOt 

where  k  can  be  chosen  to  be  independent  of  b  ^  bg.  Then, 
making  the  transformation  t  -►  bt,  we  obtain: 

^  1  k^  t‘“  e‘''‘^^|w|  2*  t  >  0 
Y“  ' 

(2.6) 

le‘^^w|^^<  k,  t  >0 

with  k|  >  kb"“,  for  w  e  and  0  <  a  <  1 ,  b  ^  bg. 

Let  us  now  return  to  the  full  equations  (1.1),  (1.2)  and 
consider  first  the  linear  case  f  *  0.  We  are  assuming  that 

|!!l-'djej{dj,x)  *  ej(x), 

for  j  =  1 ,2,. . .  ,N. 

Let  ’^i(dj)  be  the  first  eigenvalue  of  (2.1)  with  (b,-a(b,*)) 
=  (dj ,ej(dj. ,• )  and  <f^(dj,x)  be  a  normalized  eigenfunction  and 
V  >  0  such  that  the  second  eigenvalue  -  ^O' 

If  d  =  min(d^ ,. . .  ,dj^) ,  then 


Ajj  *  dl39( 


(2.7) 


?  »  diag(5i»...»5f^) 


Let 


(2.8)  = 


0-|(d^) 


...  0 


<^-|(dN) 


Xa  ^  (Y“)^ 


and  decompose  X  as 


X"-*  =  u  €“  u-*- 

a 

(2.9)  U  =  {(^  =  ♦p  v:  V  e  } 

=  {;  e  X"*:  <V,<>  =  0  for  all 


where  <.  .:>  =  f  .^(x);(x)dx. 
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If  u(t,x)  is  a  solution  of  (1.1),  (1.2)  written  as  (1.6), 
then  V,  w  satisfy  (1.7).  The  equation 

3w/3t  *  Daw 

3w/3n  =  Ew  on  3fi 

generates  an  analytic  semigroup  on  which  from  (2.6) 

satisfies  the  estimates 

|T.(t)w|  ,  <  kt'“  e'^'^^^lwl  .  .  t  >  0. 
a  U 

|TD(t)wl^  1  ke'^^^^jwl^,  t  >  0. 

a 

Under  the  hypothesis  of  the  theorem,  the  equation  (1.5)  has 
a  compact  attractor  A  As  in  [5],  there  is  a  Lipschitz  continuous 

?  ,T 

Liapunov  function  V:  such  that,  for  any  v  e  and  A  =  A  , 

y  -9  >  ^ 

(i)  V(v)  =0  if  V  e  Ag 

(ii)  a(d(v,Ag))  <  V(v)  <  b(d(v,Ag))  where  a(r)  is 

continuous,  nondecreasing,  a(r)  >  0  if  r  >  0  and  b(r) 
is  continuous,  b(0)  =  0. 


-V(v)  where 


“  TTm^^+h"^[V(v(h,VQ))  -  V(vq)] 


with  v(t,VQ)  being  the  solution  (1.5)  through  Vq  at  t  =  0. 

For  any  c  >  0,  let  =  {v  €  :  V(v)  <  c},  ®  Cz  V^. 

From  (ii),  V_  is  compact  for  any  c  >  0.  Assume  now  that  3/4  <  a 

W 

00  |J 

<  1  so  that  L  {fi,P  )  is  continuously  embedded  in  X“.  As  in 

Hale  [5  ],  we  will  show  that,  for  any  c^  <  c,  there  are  6«  n, 

n  >  6  >  0  and  bQ  such  that  Vq  €  ,  |wq |  ^  <  6,  b  ^  bQ  imply 

1  X 

that  the  solution  (v(t),w(t))  of  (1.7)  through  Vq,  Wq  satisfies 

u(t)  =  w(t)  €  W.  for  t  lO,  where 

u  c^ 


Let  u(t,‘)  be  a  solution  of  (1.1),  (1.2)  and  let  utt,*)  *  tp(*)v(t) 
+  w(t,*).  v(t)  €  p^  ,  w(t,*)  €  u^.  Then: 

dv/dt  =  P(v,w) 


(2.10) 


9w/3t  =  Diw  +  Q(v,w)  in  n 


where 


3w/3n  =  Ew  on  3f5 


P(v,w)  =  .'.qV  +  ^r)(x)f(«r)(x)v  +  w(x))dx 


Q(v,w)(x)  =  +  w(x))  -  w(y))dy- 


Transforming  v  ->■  |npv,  and  using  the  variation  of  constants 
formula,  we  rewrite  this  as 


dv/dt  =  cv  +  f(v)  +  [|fir*^( |fi|’^v,w)  -  cv  -  f(v)],  v(0)  *  v. 


(2.11) 


'  ft 

w(t)  =  Tp(t)wQ  +  Tp(t-s)Q(|fil**v,w)ds 


Since  a  >  3/4,  C  X“  implies  <*>  e  L“(fi,F^)  and  there  is  a 

constant  k  such  that  kl  <  kU. I  .  Hence,  there  is  an  n 

L*  ”* 

such  that  Vq  €  and  [Wq]  ^  <  n  imply  Vq  +  Wq(x)  e  fo 
1  X 

all  X  £  Q.  Next,  observe  that 


I  |fir''^P(  k|‘'^V,w)  -  CV  -  f(v)l  =  |(Ap  -  5)v 

+  "  ifjr^t(v+w)] 

+  [f(v+w)  -  f(v)]|  <  u-i  +  M^Iw] 

k  L 


i  *  "1 


where  =  sup{lf'(v)l:  v  €  and  =  V|(d) 


0  as  d  -*• 


In  the  same  way,  we  obtain 


|Q(  k|^v,w(x))  I  <  |f($p(x)  |nl^v+w(x))  -  f{v+w(x))j 

+  lf(v+w(x))  -  f(v)|  +  kn(x)f  $n^(tnV+w)  -  Iri'M  f(v+ 


1-1 


+  IfiT'l  |f(v+w)  -  f(v)|  <  ^2!^!*'^  *  2M  |w] 
Jj2  "  ^  L 


1  2McMw|^^  +  P2l^l 


-I5 


and 


lQ(|n|%,w)|j^  l2M^klnP|wl  ^  +  ^2 

X 


where  ug  “  ^2^^^  -*  0  as  d  ->  «. 


Then,  for  v(t),  w(t)  satisfying  (2.10),  if  v(s)  e  V  , 

*^1 


lw(s)|  <  n  for  0  ^s  ^  t,  relations  (2.11)  imply  that 


'J(v(t))  <  -V(v(t))  +  k2  I|5ir^P(  lof'^v.w)  -  i;v  -  f(v)] 


<  -V(v(t))  +  k2kMj.jwj  ^  +  k2'.;i 

X 


2(t)  ^  k-jC 


-(w-c)dt 


z(0)  +  Z\n\\k 


rt 


(t-s)’^e‘^’“'"^'^^^‘^^z(s)ds 


+  kT^2®'^'^j^(t-s)'V"‘^^^'^^ds 

<  k^e’^''’''^'*^z(0)  +  2!n|^k^kM^L(vd)'^"''y(t)  +  k^u2'-(’^d)"‘''d'*‘^ 


where  k-  is  the  Lipschitz  constant  for  V  on  V  ,  0  <  c  <  •. , 

•  OC' 

z(t)  =  lw(t)]  y(t)  =  sup{z(s),  0  <  s  <  t}  and  L  =  [  s’-'e‘^^’'^'^^cJ 


-1 


•'0 


Choosing  d  ^  dQ  such  that  f  =  1  -  2ln|'^k.|k  M^L(ud  )^‘  >  0,  we  have 


y(t)  1  k^;-^[z(0)  + 
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,a-l 


where  ^2  =  L(pd)  vj2-  implies 


(2.12a) 


If  v(s)  e  V  and  lw(s)|  <  n  for  0  <  s  ^  t,  we  also  have: 

Cl 


(2.12b)  ^(v(t))  <  -V(v(t))  +  l<2(kM^n  +  M■^). 


Let  V(v)  >  D  >  0  for  v  e  V.'vV  ,  and  choose  n  >  0  and  dr,  so 

—  C  C-j 

that: 


P  -  k2(k  M^n  +  V2)  >  0- 


Then  choosinq  6  and  d  ^  dg  so  that 


k, e  (6  +  Po)  n» 


relations  (2.12)  imoly  that,  if  v,,  €  V  ,  lWo|  <  6,  then  v(t)  e  V  , 

^  ^1  ^  '“I 

|w(t)l  <  Ti  for  all  t  >  0.  This  finishes  the  proof  of  the  claim. 

Now  let  A  denote  the  following  subset  of  U: 


A  =  {OgV  :  V  e  A.^^K 


Then  B  {Ug  =  4gVg  +  Wg,  Vg  €  ,  |Wgl  <  6}  is  a  bounded 

1  ^  .j. 

neighborhood  of  A  and  the  positive  semiorbit  (B)  is  in  W 


and  is  precompact.  Moreover,  the  estimate  (2.12a)  implies  that, 

if  d  is  sufficiently  large,  there  exists  a  neighborhood  N(a) 

such  that  Be  N(A)  and  the  u:- limit  set  u)(B)  of  B  satisfies  u(B)cV(A), 

u(B)  *  Ct  Tjj  £(t)B.  Then,  by  the  results  in  [7  ,  theorem  5.3], 

we  obtain  the  existence  of  a  comoact  attractor  Ap  in  N(A). 

We  want  to  prove  that  is  a  manifold  which  is  essentially 

the  same  as  A  r.  To  do  this,  we  construct  a  local  integral  manifold 

? » ' 

w  =  h(v,D),  V  e  V.  of  (1.7).  If  one  extends  f\V  :  V  ■*  F*'*  to 

C*j  t  t 

all  F*'*  in  such  a  way  as  to  obtain  a  uniformly  C^’^(F^,  F^)  function, 
the  referred  to  integral  manifold  must  satisfy  the  differential  integral 
equation 


(2.13) 


dv/dt  =  P(v,h(v»D)),  v(0)  =  Vq 


(2.14) 


rO 


h(vQ,D)  =  Tp(-s)Q(v(s),h(v(s),D))ds 


where  v(s)  =  v(s,VQ,h)  is  the  solution  of  (2.13)  through  Vq.  In 
a  more  or  less  standard  but  technical  way  (see,  for  example,  Henry 
[9]  for  PDE's  and  Hale  [8]  for  ODE's),  one  uses  the  contraction 
mapping  on  an  appropriate  set  of  C^’^(F^  ,U^)  functions  h(»,D)  to 
show  that  dQ  may  be  chosen  large  enough  so  that  (2.13),  (2.14)  has 
a  solution  h(v,D)  for  d  i  cIq.  Restricting  to  v  €  ,  we  obtain 

the  positively  invariant  manifold 


Mr.  =  {u  =  CnV  +  h(v,D),  V  e  V  } 
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Moreover,  h(v,D)  and  its  derivative  approach  zero  as  d  -►  «,  the 
manifold  Mq  is  exponentially  asymptotically  stable  and  the  flow 


on  it  is  given  by  (v(t ),h(v(t),D))  where 
is  a  solution  of  the  ODE  (1.9).  Since 


v(t),  v(0)  *  v«  €  V 
.u  1 

c  and  tjj  -►  , 


it  follows  that  the  vector  field  in  (1.9)  appropriately  rescaled 
(v  -*■  IfijS)  uniformly  approaches  g(v,c)  »  cv  +  f(v)  in  (1.5). 


Thus, equation  (1.9) has  an  attractor  Bp  in  .  The  attractor 
A_  f  of  (1.5)  isupper  semicontinuous  (see  [  7])  and  so  lim  Br, 
c  The  attractor  A^  for  (1.7)  is  given  by 


Ap  =  {u:  u  =  4»pV  +  h(v,D),  v  €  Bp} 


Since  h(v,D)  ■*  0  as  d  this  completes  the  proof  of  Theorem  1.1 


3.  A  scalar  example  In  one  dimension 


Consider  the  scalar  equation 

(3.1)  3u/3t  *  +  f(u),  0  <  X  <1, 

with  the  boundary  conditions 

Uy  -  ^  at  X  »  0 

(3.2) 

Uj^  +  e^u  =  0  at  X  =  1 

where  d  >  0  and  Ci  =  Ci(d)  are  constants  such  that 

J  J 

(3.3)  d  ej(d)  j  *  0,1,  as  d 

J  J 

The  equation  corresponding  to  (1.5)  is 

(3.4)  dv/dt  =-(i;o  * 

Theorem  1.1  asserts  that.  If  (3.4)  has  a  compact  attractor 
then  (3.1)  has  a  compact  attractor  if  d  >  dQ,  dg  suf¬ 
ficiently  large.  Furthermore,  if  the  flow  for  the  scalar  equation  is 
structurally  stable,  then  the  flow  on  is  structurally  stable  and 
is  equivalent  to  the  one  on  A^  The  attractors  Aj  and  A,  , 


are  close  after  A  ,  is  appropriately  embedded  inttie  function 

t  »T 


space. 

The  purpose  of  this  section  is  to  indicate  how  the  perturbed 
boundary  conditions  affect  the  flow  on  the  attractor.  To  do  this, 
let  us  impose  some  additional  conditions  on  f  which  will  ensure 
that  (3.4)  has  an  attractor  for  all  Cq,  .  More  precisely,  suppose 
that  there  is  a  Vq  >  0  such  that 


(3.5) 


vf(v)  <  0  for  Ivl  >  \ff 


f(v)/v  as  1  V|  •+  « 


For  any  Equation  (3.4)  has  a  compact  attractor  ^  which 
is  an  interval  K  ^  =  [a,B]  where  a,  6  are  the  extreme  zeros  of 
(^0  +  ^(v). 

As  a  particular  illustration  suppose  f(v)  has  five  simple  zeros 
as  shown  in  Figure  1. 


Figure  1 


For  Neumann  boundary  conditions  and  d  i  dQ,  the  attractor  Aq  ^ 

Is  the  segment  [ai*ag]  and  the  flow  on  the  attractor  is  the  one 
indicated  in  Figure  1.  By  varying  c,  it  is  clear  that  one  can 
obtain  each  of  the  attractors  shown  in  Figure  2.  The  perturbatiesns 
in  the  boundary  conditions  induce  large  changes  in  the  dynamical 
properties  of  the  system. 


Figure  2 

The  flow  on  the  attractor  in  the  function  space  for  (3.1),  (3.2) 
is  the  same  as  the  ones  indicated.  However,  the  equilibrium  points 
are  not  constant  functions.  The  exact  equilibrium  points  are  not  easy 
to  compute  since  one  must  solve  a  nonlinear  boundary  value  problem. 
Figure  3  in  the  (u,u^)-plane  Indicates  the  approximate  location  of 
these  equilibria  for  the  case  eq  >  0,  e-j  >  0.  The  line  Lq  =  {(u,Uj^) 
ux  -  EnU  =  0}  and  the  line  L,  *  {(u,u„):  u„  +  e,u  =  0}. 
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4.  Two  examples  In  genetic  repression. 
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and  transfer  through  the  membrane  bounding  u.  This  model  was  dis¬ 
cussed  in  the  one  dimensional  case  by  Mahaffy  and  Pao  [11]  and  in  the 
general  case  by  Busenberg  and  Mahaffy  [1]. 

Let  u^  and  v^,  i  *  1,2  be  respectively  the  concentrations  of 
the  mRNA  and  repressor  protein  in  compartments  u  and  (Tno.  The 
nucleus  u  is  considered  as  a  well  mixed  compartment  containing  mRNA 
whose  concentration  u^  is  transcribed  from  the  gene  at  a  rate  depending 
on  the  repressor  protein  v^.  The  mRNA  leaves  u  and  enters  the  cyto¬ 
plasm  fixj  where  it  diffuses  and  interacts  with  the  ribosomes.  Through 
the  delayed  process  of  translation,  a  sequence  of  enzymes  is  produced 
which  in  turn  produces  a  repressor  V2.  This  end  product  diffuses  back 
to  u  where  it  inhibits  the  production  of  u-j.  The  cell  wall  is  a  bar¬ 
rier  through  which  neither  biochemical  substance  can  pass  (see  Figure  4). 

Let  w^  *  col(u.|  ,v-j) ,  W2  =  col(u2.V2)»  B  »  diag  (bi,b2),  A  ® 
diag(a^,a2),  D  =  diag(d^,d2),  C  =  ^2^  =  Cq,  c^j  »  0  otherwise, 

S  =  diag(6p62)«  =  col(f(v^),0),  where  b^,  a^,  dj,  6^,  Cq  are 

positive  constants  and  f(v)  is  a  decreasing  function  of  v  of  the 
form  [l+v*^]"^  where  p  >  2  is  a  constant  and  v  €  F.  Also,  let 
r^  10,  r2  lO,  be  given  constants.  With  this  notation  the  two  com¬ 
partment  model  is 

dw,(t)  r 

—  =  F(w^(t-r^))  -  Bw^(t)  A  [W2(t,y)  -  w^(t)]dy 

3w2(t,x) 

— - =  DAw2(t,x)  -  Bw2(t,x)  +  Cw2(t-r2,x),  x  €  rrv, 


(4.1) 
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with  the  boundary  conditions 


(4.2) 


3W2(t,x) 

— gfj —  =  -eCwgCt.x)  -  w.|{t)]  on  au) 


(4.3) 


3W2(t,x) 


an 


0  on 


an. 


where  n  denotes  the  outward  normal  to  ar  with  r  *  jtsu. 

The  flow  defined  by  these  equations  leaves  the  region  = 

{Uj  >  O.Vj  >  0,  j  =  1,2}  positively  Invariant  (see  Mahaffy  and  Pao 

[11]). 

To  apply  the  methods  In  the  proof  of  Theorem  1.1,  It  Is  convenient 
to  make  a  transformation  of  variables  In  (4.1),  (4.2),  (4.3): 


w^(t)  -►  W^(t)  W2(t,x) 


W2(t,x)  +  W^(t). 


The  new  equations  are 


<lWi(t)  , 

(4.4)  _g!_  =  F(W  (t-r^))  -  BW,(t)  +  A  W2(t,y)dy 

*  ** 

dW  { t  X ) 

"ir  '  D4W2(t,x)  -  BW2(t.x)  -  F(W,(t-r^)) 


+  CW2(t-r2,x)  +  CW^(t-r^,x)  -  A  W2(t,y)dy 

•  rt. 


with  the  boundary  conditions 


We  also  need  the  space  in  which  we  will  consider  (4.4),  (4.5),  (4.6). 

p  p 

Let  X  «  )  and  x“,  0  ^  a  <  1  be  the  usual  fractional  power 

space  discussed  before  and  associated  with  the  Laplacian  with  boundary 
conditions  3u/3n  «  -bu  on  3u,  au/3n  =  0  on  an.  Let  X^,  x“  be 
the  restrictions  to  the  positive  functions.  Equations  (4.4),  (4.5), 

(4.6)  then  define  a  C^-semigroup  on  C([-r,0],Rj)  x  C([-r,0],Xj). 

As  remarked  earlier,  the  method  of  proof  of  Theorem  1.1  applies 

to  this  set  of  equations  when  d  *  min(di,d2)  is  large  and  the  Bj 
are  functions  of  d^,  so  that 

(4.7)  l3a)|djBj  -►  Cj  as  dj  -►  » 

The  resulting  set  of  ordinary  differential  difference  equations  which 
are  the  analogue  of  equation  (1.5)  are 

dW,(t) 

- F(W^(t-r^))  -  BW^(t)  +  A^W2(t) 

(4.8) 

dW2(t) 

-  -cW2(t)  -  BW2(t)  +  CW2(t-r2) 

+  CW^(t-r2)  -  A^W2(t)  -  F(W^(t-r^) 
where  :  »  diag(  .j.-p),  =  A  , 
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Mahaffy  and  Pao  [11]  have  shown  that,  for  certain  values  of  the 
parameters  in  (4.8),  there  is  a  globally  attracting  stable  equilibrium 
point.  Therefore,  the  attractor  is  a  single  point.  They  show  also 
that  a  stable  generic  Hopf  bifurcation  occurs  at  this  equilibrium  point 
with  the  delays  being  used  as  the  bifurcation  parameter.  Busenberg  and 
Mahaffy  [1]  show  that  (4.4),  (4.5),  (4.6),  for  the  diffusion  coefficients 
large  enough,  must  also  have  a  generic  Hopf  bifurcation  at  approximately 
the  same  value  of  the  delay  parameters.  This  is  proved  by  comparing  the 
characteristic  equations  of  the  linear  variational  equations  of  (4.8) 
and  (4.4),  (4.5),  (4.6)  near  their  respective  equilibria. 

Let  us  now  discuss  the  implications  of  our  theory  to  this  situation. 
Suppose  c  is  fixed  and  that  (4.8)  has  a  compact  attractor  o/  ^  ^ 

;,ri,r2 

in  C([-r,0],F^).  Under  this  assumption  and,  if  (4.7)  is  satisfied, 
one  can  conclude  that  (4.4),  (4.5),  (4.6)  has  a  compact  attractor 
^  r  r  r  C([-r,0],x“)  if  3/4  <  a  <  1  and  d  >  d«  where  d^ 

u»4»ri,r2  ♦  -  u  u 

is  sufficiently  large.  Furthermore,  the  flows  on  the  two  attractors 
are  equivalent  if  the  one  on  is  structurally  stable.  Moreover, 

from  the  remarks  made  in  the  introduction,  the  Hopf  bifurcation  with 
respect  to  the  delays  referred  to  above  occurs  onja/L  ^  ^  if  d  is 
sufficiently  large. 

The  following  fact  is  also  true  even  though  it  will  neither  be  used 
here  nor  proved  in  detail.  Since  the  only  nonlinearity  F(v)  in  the 
equation  is  bounded  together  with  its  first  and  second  derivatives,  one 
can  choose  a  global  Liapunov  function  V:C([-r,0] ,F^)  -*■  F  for  the 


as  a  quadratic  form  outside  a  large  ba11.  This 


attractor  j}/  ^  ^ 

^  *  £ 

permits  one  to  obtain  a  global  attractor  for  (4.4),  (4.5),  (4.6). 

The  second  example  that  we  are  going  to  consider  Is  again  a  model 
In  genetic  repression  with  the  same  features  as  the  one  above,  except 
that  now  the  process  of  production  of  the  repressor  takes  place  In  the 
membrane  bounding  the  cell,  looked  upon  as  a  third  compartment  S.  Thi 
model  Is  also  discussed  In  the  one  dimensional  case  by  Busenberg  and 
Mahaffy  [1].  We  need  some  new  notation.  Let  u^,v^,  1  =  1,2,3  denote 
the  concentrations  of  mRNA  and  repressor  respectively  in  the  compart¬ 
ments  w,  and  the  outside  membrane  si.  The  compartments  w  and 
si  are  considered  well  mixed.  Let  w^  =  col(u^,v^),  1  =  1,2,3,  and 
B,  A,  D,  C,  e,  F  be  the  same  as  above.  Let  also  'A  -  diag(^|,a2), 

?  *  diag(6^,F2)  with  aTj,  positive  constants.  Then  the  three  com¬ 
partment  model  is 


(t) 

dt 


=  F(w^(t-r^))  -  Bw^(t)  +  A 


[w2(t,x)  -  w,(t)]dx 


3W2(t,x) 

(4.9)  — - =  DAW2(t,x)  -  B  W2(t,x),  x  €  Slsu 


dw3(t) 

dt 


=  -B  W2(t)  +  C  W2(t-r2,x)  +  A 


[W2(t,x)  -  W3(t)]dx 


with  the  boundary  conditions 


3w2(t,x) 

— ^ - =  -6[w2(t,x)  -  w^(t)]  on  3;.. 

(4.10) 

9W2(t,x) 

— ^ - =  -6twp(t,x)  -  w^(t)]  on  30 


The  flow  defined  by  these  equations  leaves  the  region  F”  = 
{u^  >  0,  >  0,  i  =  1,2,3}  positively  invariant.  Again,  it  is 

convenient  to  make  a  transformation  of  variables  in  (4.9),  (4.10) 
in  order  to  obtain  homogeneous  boundary  conditions.  Let  H  = 
col (hph2,R^ ,^2)  be  the  solution  of  the  following  boundary  value 
problem  with  nonhomogeneous  boundary  conditions: 

aH  =  0,  in  r  =  ns*) 


+  6H  =  Y,  on  sr  =  an  u  a. 


where  6  =  diag(6^  .Bg.fei  .62) »  v  =  co1(Bi  ,62*0,0)  ^  ~ 

diag(6.|  ,$2*61  ,62)s  7  =  001(0,0,61,62)  O''  Then,  let  H  =  diag(hi,h2), 
FT  =  diag(hi,R’2)  and  introduce  in  (4.9),  (4.10)  the  change  of  variables 


Wi(t)  Wi(t),  W2(t,x)  W2(t,x)  +  H(x)Wi(t)  +  H(x)W3(t),  W3(t)  1^3(1) 

If  we  let  6  =  f  H,  Gp  =  f  H,  G",  =  f  H ,  G.  =  f  H,  A,  =  la^lA 
_  ,  Jau)  1  '  ' 

and  Ai  =  1 3^2 [A,  the  new  equations  are 

dW,  r 

=  F(Wi(t-ri))-(3+Ai-AG^)Wi(t)  +  A  W2(t,y)dy  +  AG^W3(t) 

9W2(t,x) 

at  =  D:.W2(t,x)  -  B  W2(t,x) 


(4.11) 


dW3(t) 


-  H(x){F(Wi(t-ri))-(Ai-AG^)Wi(t)+Aj  W2(t.y)dy  +  AG^W3(t)} 

-  H(x){CW3(t-r2)-(Ai-AG^)W3(t)+A[  W2(t,y)dy  +  AG^Wi{t):- 


=  C  W3(t-r2)-(B+Ai-AG„)W3(t)+Aj_^W2(t.y)dy  +  AG,Wi(t) 


' II  •  *  "'v  •  V «  L iiu  lij  • 
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with  the  boundary  conditions 


(4.12) 


aWpCt.x) 

- ^ — =  -B  W2(t,x),  on  8co 


3W2(t,x)  _ 

- — =  -S  W2(t,x),  on  3fi. 


✓  4 

These  equations  define  a  C  semigroup  on  C{[-r,0],P^)  x  x“  and 
we  can  apply  again  the  method  of  proof  of  Theorem  1.1  when  d  = 
min(d^,d2)  is  large  and  are  functions  of  so  that 


(4.13)  |r|'''dj(l3.,l6j  +  |an|6j)  -  cj 


as  d 


If  we  denote  by  ?  =  diag(  ?-j  ,^2)»  G  =  ir|"M  H  and  IT* 

J  r 


'  H, 


the  resulting  set  of  ordinary  differential  difference  equations  are 


dW^(t) 


^  =  F(W^(t-r^))  -  (B+A^-AGjW^(t)  +  A^W2(t)  +  AG^W3(t) 


dW,(t) 


=  -(t  +  B  +  GA^  +  GA^)W2(t) 


(4.14)  -  GF(W^(t-r^))  +  [GA^  -  GAG^  -  GAGjj]W^(t) 


-  GCW3(t-r2)  +  -  GAGj^  -  GAG^,]W3(t) 


dW,(t) 


— =  CW3(t-r2)  -  (B+A^-AGj-pW3(t)  +  AW2(t)  +  AGj^W^(t) 


This  set  of  equations  define  a  C^-semigroup  on  C([-r,0],F°)  and 
it  is  this  flow  that  should  be  compared  to  the  one  defined  by  the 
equations  (4.9),  (4.10). 

As  remarked  earlier,  Busenberg  and  Mahaffy  [1]  considered  (4.9), 
(4.10)  in  one  space  variable  and  made  a  comparison  to  a  set  of  equations 
which  they  claimed  would  correspond  to  a  well  mixed  model.  These  equa¬ 
tions  were  analogous  to  but  not  the  same  as  (4.14). 

For  equations  (4.9)  and  (4.10)  to  correspond  to  the  physical  problem 
there  must  be  some  conservation  laws  and  this  imposes  relationships  be¬ 
tween  the  constants  A,  6,  ?.  Once  these  conditions  are  imposed, 
equations  (4.14)  in  the  one  dimensional  case  take  a  simpler  form.  For 
special  values  of  the  parameters,  one  obtains  the  equations  of  Busenberg 
and  Mahaffy  [1].  The  higher  dimensional  case  always  contains  terms  that 
depend  on  the  shape  of  the  region  in  a  significant  way.  It  is  obvious 
that  these  equations  need  to  be  investigated  in  more  detail. 
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